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ABSTRACT 


* 

The  modal  method  developed  by  the  authors^  applicable  to  the 
dynamic  analysis  of  supported  noncircular  cylindrical  shells  is  extended 
in  the  present  work  to  deal  with  (but  not  restricted  to)  the  free  vibra- 
tion problem  of  unsupported  noncircular  cylinders.  The  required  modifica- 
tion includes  the  enforcement  of  the  edge  conditions  at  the  free  ends  of 
a finite  cylindrical  shell  which  are  unsatisfied  by  the  modal  functions 
owing  to  the  presence  of  the  variable  curvature  terms.  Such  conditions 
are  posed  as  additional  constraints  by  way  of  the  well-known  formalism 
of  Lagrange  multipliers.  In  addition,  it  was  found  that  the  inclusion 
of  the  lowest  modes,  commonly  approximated  by  the  Rayl e igh-Love  modes, 
is  essential  to  the  completeness  of  the  eigen-function  representation. 

The  validity  of  the  proposed  procedure  of  analysis  is  illustrated 
by  its  application  to  the  solution  of  the  free  vibration  problem  of  oval 
cylindrical  shells  with  free  ends. 


L — 


Introd uct ion 


In  an  earlier  article*,  the  present  authors  presented  an 

▼ 

analysis  of  the  natural  vibration  of  noncircular  cylincrical  shells. 

In  addition  to  the  discussion  on  the  contrast  (as  well  as  the  similarity) 
of  dynamic  characteristics  between  circular  and  noncircular  cylindrical 
shells,  emphasis  was  placed  on  the  development  of  a modal  method  suitable 
for  application  to  problems  in  dynamics  of  noncircular,  finite,  cylindrical 
shells  having  various  types  of  simple  and  clamped  supports. 

, The  most  profound  difficulties  in  analysis  involving  cylindrical 

shells  with  variable  radius  of  curvature  stem  from  the  fact  that  the  depend- 
ency of  deformation  (and  velocity)  of  the  shell  element  upon  the  two  in- 
plane coordinates  is,  in  general,  not  separable,  as  is  the  case  of  cylind- 
rical shells  of  constant  radius.  Such  complexity,  compounded  w i th  the 
inherent  difficulties  accompanying  the  rigorous  enforcement  of  the  boundary 
conditions  corresponding  to  various  types  of  end  conditions  of  a finite 
cylinder,  nevertheless,  can  be  handled  by  the  aforementioned  modal  approach, 
as  shown  in  Ref.  1.  Such  a procedure  involves  expanding  the  middle  surface 
displacement  parameters  in  series  form  on  a functional  basis  constituted 
by  the  modal  functions  derived  from  the  exact  solution  of  the  correspond- 
ing free  vibration  problem  of  a circular  cylindrical  shell.  The  validity 
of  this  procedure  is  justified  by  the  fact  that,  within  the  framework  of 
Donnel-type  shell  theory  for  the  class  of  problems  considered  in  Ref.  1, 
the  boundary  conditions  prescribed  for  the  oval  shell  are  identically  sat- 
isfied, inasmuch  as  the  modal  functions  sat-isfy  the  required  conditions  of 
the  circular  cylinder.  The  same  conclusion  is  also  valid  for  most  cases 
when  Sanders'  shell  theory  is  employed.  Even  the  exceptions  may  be  over- 
looked, since  the  error  induced  proves  to  be  negligibly  small.  Thus,  this 

. 
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2. 

procedure  permits  the  concentration  of  effort  to  cope  with  the  problem 
of  circumferential  modal  coupling. 

This  method  of  analysis,  developed  as  a variation  of  the  modal 
method,  must  still  be  modified  further  to  be  used  in  the  dynamic  analysis 
of  oval  cylindrical  shells  with  free  ends.  In  the  present  work,  as  a 
sequel  to  Ref.  1,  such  modifications,  together  with  the  characteristic 
behavior  of  the  free-free  oval  cylindrical  shell  resulting  from  the  applica- 
tion of  the  modified  method,  is  discussed  in  detail.  Primarily,  the  mod- 
ification involves  the  posing  of  the  unsatisfied  boundary  conditions  as  an 
additional  restraint  on  the  variational  problem  within  the  framework  of 
the  modal  expansion. 


The  Energy  Functional 

In  Ref.  1,  a variational  equation  of  motion  based  on  Hamilton's 
principle  was  obtained  in  the  form  of 


6F  = 0 


0) 


in  which 


6F  = f f {[u  +v(v  +w/r)]6u  + (l/2)(1-v)(v  +u  )(6v  +6u  ) 
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In  this  expression,  the  middle  surface  displacement  parameters  u,v,  and  w, 
the  coordinates  x and  s,  and  the  radius  of  curvature  r = r(s)  are 
all  dimensionless  quantities  normalized  with  respect  to  the  shell  thick- 
ness h and  an  average  radius  rQ  in  the  following  manner  (see  Fig.  I): 

(u,  v,  w)  = (u",  v",  w )/h 

(x,  s,  r)  - (x",  s“,  r )/rQ  (2) 

where  u x etc.  denote  the  correspond ing  physical  parameters 

of  the  shell;  rQ  is  defined  to  be  (L^/217),  where  LQ  denotes  the  cir- 
cumferential length  of  the  cross-sectional  contour.  In  addition,  v is 
Poisson's  ratio,  is  a tracing  parameter  which  is  equal  to  1 for  Sanders' 

shell  theory  and  equal  to  zero  for  Donnell's  theory.  The  constant  k and 
the  parameter  X are  defined  as 

k = (1/12)  (h/rQ)2  , X = n2  = pu>-  rQ2(l-v2)/E  (3) 

in  which  x denotes  the  natural  frequency  of  the  oval  shell,  and  p and 
E are  the  mass  density  and  modulus  of  elasticity,  respectively.  The 
development  of  the  preceding  variational  equation  and  the  required  modal 
functions  were  discussed  in  Ref.  1 and  will  not  be  repeated  here. 

It  is  worth  noting  that  the  variational  functional  will  be  main- 
tained in  the  present  form,  which  offers  the  optimal  degree  of  operational 
symmetry.  However,  the  usual  integrations  by  parts,  leading  to  the  ex- 
plicit determination  of  the  equations  of  motion  and  the  natural  boundary 
conditions  is  needed  to  facilitate  the  discussion  of  the  latter.  To  this 

end,  the  quantities  5u,  6v,  6w,  and  6w  are  regarded  as  the  chosen 

» x 

virtual  displacements.  It  follows  that,  if  the  ends  of  the  oval  cylinder 
are  free  from  traction  (which  is  the  primary  concern  of  the  present  work), 


T 


f 

4. 


the  following  quantities  are  required  to  vanish  at  both  ends  of  the  shell; 
i .e . , 

u + v(v  +w/r)  = 0 * 

»x  »s 


v y + u - (3kk  /4)  (1/r)  (4w  - 3k  v Jr  + k u /r)  = 0 

5 »Xb  S >X  S ,S 


w + v (w  - k v/r ) = 0 

,xx  , s s ,s 


(4) 


w XXX  + v(w  s - k v/r)  + (1/2)  (1-v)  (4w  -3k  v Jr  + k u Jr)  c=0 

i AAA  i p 5 ) X5  >X5  S | X S y S 9S 


The  associated  equations  of  motion,  which  represent  a form  of  the  well- 
known  Sanders'  equations,  are  not  shown  here. 

Equations  (4),  in  essence,  state  that  the  axial  resultant  force, 
the  in-plane  shear,  the  bending  moment,  and  the  effective  transverse 
shearing  force,  respectively,  must  vanish  at  the  boundary.  Alternative 
boundary  conditions  for  the  various  types  of  end  supports  (simple  and 
clamped)  were  discussed  in  Ref.  1. 


Modal  Expansion  and  Kinematic  Constraints 


Similar  to  the  cases  of  the  supported  cylinders  analyzed  in 
Ref.  1,  the  displacement  parameters  u,  v,  and  w are  assumed  to  be 

N 

u = .£.  ar.x.U  (x)cos(ns) 
i=l  i in 

N 

v = .£.  y . V (x) s i n (ns ) 
i = l i n 

N 

w = X]  a.z.Wn(x)cos(ns)  (5) 


L. 


5. 


where  the  subscripted  indices  i and  n are  related  by 
, (n+l)/2  if  n is  odd 

i = j 

"(n/2)+1  if  n is  even  T (6) 

and  cr.  is  equal  to  1//2  if  n = 0 and  is  otherwise  equal  to  1. 

It  should  be  noted  that  it  has  been  assumed  that  u,  v,  and  w are 
sinusoidal  in  time  with  circular  frequency  ou  ; i.e.,  u = E ar.x.U  (x) 
cos  (ns ) s i n (cut ) , . . . , etc.  The  time  factor  sin(out)  is  deleted,  and  its 
presence  is  reflected  by  the  introduction  of  the  parameter  X , 
defined  by  Eq . (3 ) » which  appears  in  Eq.(l).  In  Eqs.(5)  the  deformations 
(as  well  as  the  corresponding  velocities)  are  assumed  to  be  symmetric 
with  respect  to  the  origin  of  the  coordinate  s ; i.e.,  s = 0.  For  the 
sake  of  brevity  in  identification,  let  the  term  "symmetric  modes"  be 
used  in  the  ensuing  discussion  to  describe  the  modes  of  deformation  asso- 
ciated with  the  expansion  given  in  Eqs.(5).  Similarly,  the  term  "anti- 
symmetric modes"  is  used  for  modes  of  deformation  given  by  similar  expan- 
sions by  replacing  the  functions  sin(ns)  by  - cos(ns),  and  by  replacing 
cos(ns)  by  sin(ns)  in  Eqs.(5).  Such  a distinction,  of  course,  is  trivial 
for  the  circular  cylinders.  Parallel  to  these  designations,  modes  that 
are  symmetric  or  antisymmetric  with  respect  to  the  mid-span  of  the  shell 
will  be  called  "longitudinally  symmetric"  or  "longitudinally  antisymmetric", 
respect ively . 

In  the  present  work,  as  in  Ref.  1,  the  nond imens ional  radius  of 
curvature  is  chosen  to  be 

r = 1/(1  + Z cos  2s)  (7) 

where  5 is  the  eccentricity  parameter,  the  magnitude  of  which  lies 

between  zero  and  unity.  Moreover,  the  quantities  U , V , and  W , 

n n n 


6. 


introduced  in  the  format  of  Eqs.(5)  (for  any  value  of  n),  represent 
the  exact  solution  to  the  free  vibration  problem  of  a circular  cylindrical 
shell  posed  by  the  variational  problem  of  Eq.(l)  for  r = 1.  These  pre- 
determined modal  functions,  by  definition,  satisfy  the  boundary  conditions 
prescribed  by  Eqs.(4)  for  r = 1.  The  existence  of  such  a set  of  exact 
solutions,  of  course,  is  well  known.  (See,  for  example,  Ref.  2)  .Detailed 

discussions  on  the  computational  aspects  of  such  a problem  may  be  found  in 

3 4 

the  works  of  Vronnay  and  Smith  , and  of  Warburton  , as  well  as  in  Ref.  1. 

Upon  substitution  of  Eqs.(5)  (or  their  "antisymmetric"  counterparts) 
into  Eq.(l),  it  follows  that,  if  the  undetermined  parameters  x.,  y.,  and 
z.  are  permitted  to  vary  independently,  a set  of  3N  algebraic  equations 
can  be  obtained,  the  details  of  which  may  be  found  in  Ref.  1.  The  presence 
of  r(s)  as  a function  of  the  coordinate  s given  in  Eq.(7)  results  in 
the  coupling  of  these  algebraic  equations.  Nevertheless,  it  rs  also  clear 
that  the  "symmetric"  and  "antisymmetric"  modes  are  not  related,  owing  to 
the  double  symmetry  possessed  by  r(s)  as  given  in  Eq.(7).  Within  each  of 
these  two  groups,  the  odd  modes  and  the  even  modes  are  uncoupled  as  well. 

The  circular  cylinder  problem  can  be  treated  as  a degenerate  case  of  this 
formulation.  In  such  a case,  coupling  exists  only  among  x.,  y.,  and  z. 
for  a given  value  of  the  subscript  n,  the  circumferential  mode  number. 

As  regards  the  boundary  conditions,  when  Eqs.(5)  (or  their  anti- 
symmetric counterparts)  are  substituted  into  Eqs.(4),  the  presence  of  r 
as  a function  of  s prevents  the  resulting  algebraic  form  of  the  boundary 
conditions  from  being  identically  satisfied. 

In  this  connection  it  may  be  advantageous  to  first  discuss  the 
simpler  case  corresponding  to  kg  = 0;  i .e . , the  "Donnell"  version  of  the 
equations  involved,  since  mast  of  the  terms  containing  the  variade  radius 


of  curvature  vanish  upon  setting  kg  = 0.  In  so  doing,  Eqs.(4)  reduce  to 


u + v (v  + w/r)  = 0 


v + u = 0 
>X  ,s 


7. 


w + vw  = 0 
, XX  , s s 


w + (2-v)  w = 0 

,xxx  ,xss 


(8) 


for  x = x , say,  where  x is  the  axial  coordinate  of  a free  end. 
e e 

Obviously,  the  third  and  fourth  conditions  of  Eqs . (8)  are  identically 

satisfied  by  the  deflection  w in  terms  of  W (x) , since  each  of  the 

n 

functions  Wn  satisfied  such  conditions.  The  second  condition  in  Eqs . (8) 
is  satisfied  if  x.  and  y.  are  related  by 


nil  (x  ) x . - V ' (x  )y . = 0 (9) 

for  all  values  of  i . In  Eq . (9)  a primed  superscript  denotes  differentia- 
tion with  respect  to  x.  Finally,  when  the  modal  expansions  are  substituted 
into  the  first  of  Eqs . (8) , it  follows  that  x.,  y.,  and  z.  must  be  re- 
lated by 


ECx.U^x  ) + vy.nV  (x  )]cos  ns  + E vz.W  (x  ) [ cos  ns  + (§/2) cos  (n  + 2) s 
ine  in  e . ine  ' ' 

I i 


+ (E/2)cos (n-2)s]  = 0 


00) 


for  the  symmetric  modes.  A similar  condition  is  arrived  at  by  employing 
the  antisymmetric  modes,  and  the  result  is  identical  to  Eq . (10)  with  the 
cosine  functions  replaced  by  the  corresponding  sine  functions.  This 
condition  can  be  satisfied  provided 


x.U^x  ) + v{y.nV  (x  ) + z.W  (x  ) + (5/2 )[  z . ,W  _ (x  ) + z.  ,W  -(x  )]}  = 0 
i n e 1 7 i nv  e i n e i+l  n+2 v e'  i-l  n-2  e J 
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for  i > 2 . 


Individual  care  must  be  administered  for  the  irregular  terms 


, _ 
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i 

when  i <_ 2 ; they  can  be  derived  from  Eq.(10).  For  instance,  the  collection 
of  coefficients  of  cos(s)  in  Eq  .(10)  yields  the  following  constraint  condi- 
tion for  i = I of  the  odd,  symmetric  group:  T 

X1UJ  (xe)  + vW,V,  (xe)  + z,wi  (xe)  + (E/2)[z1W]  (xg)  + z2W3  (xe)]}  = 0 

The  constant  a.  has  been  deleted  from  Eqs.(10)  and  (11),  but  it  must  be 
restored  to  several  terms  whenever  needed. 

Hence,  the  variational  problem  posed  by  Eq.(l)  must  be  solved  with 
the  kinematic  constraint  conditions  represented  by  Eqs.(9)  and  (11),  for  all 
values  of  index  i.  Suppose  that  the  modal  expansions  are  truncated  to  retain 
the  first  N terms  in  each  of  the  three  series  expansions;  the  unconstrained 
system  of  equations  consists  of  3N  equations  for  the  3N  undetermined  para- 
meters. (The  mode  Vg  in  the  symmetric  case  and  Ug  and  Wg  in  the  anti- 
symmetric case  may  be  retained  to  preserve  formal  symmetry  but  they  are  not 
coupled  with  the  other  modes).  Equations  (9)  and  (11)  can  then  be  utilized 
(for  instance,  with  use  of  Lagrange  multipliers)  to  reduce  the  system  of  equa- 
tions to  one  that  consists  of  N equations  for  the  N undetermined  para- 
meters z..  While  such  reduction  is  a welcome  relief  computat ionwise,  the 
same  reasoning  indicates  that  the  use  of  Sanders'  equations  in  the  analysis 
(and,  perhaps,  other  shell  theories  as  well)  would  lead  to  over  restraint, 
since,  upon  retaining  terms  in  Eqs . (4)  that  were  dropped  under  kg , the  bend- 
ing moment  and  the  effective  shearing  force  conditions  produce  two  additional 
sets  of  kinematic  constraints.  This  situation  thus  dictates  the  limitation 
of  the  present  method  of  analysis.  Nevertheless,  based  on  the  numerical 
verification  in  Ref.  1 that  a similar  degree  of  violation  of  the  edge  condi- 
tions when  Sanders'  equations  are  used  is  only  a minor  one,  it  is  believed 
that  the  use  of  Sanders'  equations  in  the  present  problem  is  also  valid  as 
long  as  the  first  two  conditions  of  Eqs. (4)  are  enforced  in  the  preceding 

A 


manner . 


(Terms  proportional  to  k in  the  second  of  Eqs.(4)  are  negligible.) 


The  Modal  Functions 

The  modal  functions  U , V , and  W , as  mentioned  in  the  fore- 

n n n 

going  sections,  are  the  exact  normal  mode  functions  for  the  free  vibration 
of  the  circular  cylindrical  shell.  It  can  be  shown  that  these  functions, 
with  subscripts  deleted,  must  satisfy  the  following  differential  equations 
(Kj=  0): 

U"  - l/2(1-v)n2U  + XU  + l/2(l+v)nV'  + vw'  = 0 

J/2(l+v)nU'+(n2  - X)V  - l/2(l-v)v"  + nW  = 0 

vU'  + nV  + k (W 1 V - 2n2Wff  + n4W)  + (1  - X)W  = 0 (12) 

and  the  following  boundary  conditions; 

U'  4 v(nV  + W)  = 0 , V'  - nU  = 0 , W"  - vn2W  = 0 

W'"  - (2-v)n2W'  = 0 for  x = +L/(2ro)  (13) 

where  L denotes  the  length  of  the  shell. 

To  solve  the  above  boundary  value  problem,  the  mode  functions 
U,  V,  and  W are  assumed  to  have  the  exponential  form 

8 

(U,  V,  W)  = Re  k|j  (ak,  bk>  cR)  exp(Bkx)  (14) 

The  procedure  to  determine  the  amplitude  ratios  ak/ck>  and  bk/ck  , 

2 

and  the  axial  wave  numbers  0k  has  been  extensively  discussed  . In 
particular,  a computational  procedure  proposed  by  Vronnay  and  Smith^  is 
judged  most  convenient  and  is  adopted  in  the  present  work.  The  natural 
frequencies  and  the  associated  characteristics  can  thus  be  determined 


10. 
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n ume  r i ca  I 1 y . 

It  is  interesting  to  note  that  the  issue  of  the  lowest  modes 
of  free  vibration  is  somewhat  obscure  in  the  literature.  The  so-called 
Rayleigh  and  Love  modes,  labeled  as  inextens iona 1 , are  generally  regarded 
as  approximations,  inasmuch  as  they  satisfy  neither  the  equations  of  motion 
nor  the  boundary  conditions.  As  a consequence,  the  modes  with  frequencies 
higher  than  that  of  the  "inextens iona I modes"  are  recognized  as  the  first 
modes . 

However,  if  the  strain-displacement  equations  of  the  Sanders' 
shell  theory  are  integrated  for  zero  strains  and  zero  edge  forces,  two,  and 
only  two,  sets  of  exact  solutions  do  exist,  both  occurring  when  the  circun- 
ferential  wave  number  n is  equal  to  1.  They  are 


and 


Vj  C^x  , Wj  - Cj.x 


05) 


in  which  C's  are  constants.  These  expressions,  being  precisely  the 
Rayleigh  and  Love  modes,  respectively,  for  circular  cylindrical  shells, actu- 
ally represent  rigid  body  translation  and  rotation.  Although  these  modes  are 
trivial  solutions  for  n = 1 , they  are  exact,  nonetheless.  This  observa- 
tion, together  with  the  experimental  evidence  obtained  by  Sewall  and 
Naumann^,  indicates  that  the  "Ray  1 e i gh-Love"  modes  do  exist,  although  they 
might  be  different  from  their  original  approximate  version  given  by  Eq  s . (15) 
and  they  might  not  be  inextens iona 1 in  nature.  In  fact,  the  attempt  herein 
to  determine  these  modes  numerically  is  successful. 

The  characteristic  wave  numbers  0^  , determined  from  the  secular 

equations  for  the  amplitude  ratios  a^/c^  anc*  f°r  a 9'ven  value  of 

2 

X , may  be  categorized  in  the  following  fashion.  The  four  roots  for  6^ 


A 


can  be  arranged  in  numerical  order  according  to  their  absolute  values. 

2 2 2 2 

As  X increases  from  zero, 8j  , 8^  » 8^  ,8^  form  two  complex  conjugate 

2 2 

pairs  (hereby  named  zone  1).  As  X increases  further,  8j  and  8^ 

2 2 

are  real  and  positive,  while  8^  and  remain  complex  conjugates 

2 

of  each  other  (zone  2).  As  X continues  to  increase,  8j  becomes  real 

2 2 2 
and  negative  while  is  real  and  positive;  and  8^  and  8^  are 

complex  conjugates  (zone  3).  When  the  analysis  is  based  on  Donnell's  shell 

theory,  the  afore-mentioned  secular  equation  is  quite  simple  and  such  zonal 

behavior  can  be  established  without  much  difficulty.  For  the  shell  model 

selected  for  numerical  study  in  the  present  work,  Rayleigh  modes (long i tud in- 

ally  symmetric)  are  found  in  zone  2,  while  Love  modes  (longitudinally 

antisymmetric)  lie  in  zones  2 and  3.  In  any  case,  they  are  not  rigid  body 

motions,  and  they  are  exact.  As  X increases  beyond  these  modes,  the  modes 

commonly  labeled  as  m = 1,2,  . . . (m  being  the  nominal  axial  wave  number) 

are  found  in  zone  3. 

While  Eqs.(15)  are  not  exact  for  modes  other  than  n = 1,  they  can 
be  used  to  determine  the  approximate  eigenvalues  of  the  Rayleigh-Love  modes 
to  provide  some  reference  for  the  numerical  computation.  Substitution  of 
the  mode  shape  functions  given  in  Eqs.(l5)  into  Eqs.(5)  in  place  of  those  of 
the  nth  mode  and  the  subsequent  introduction  of  the  resulting  expressions 
into  Eq.(l)  with  r = 1 yields  the  approximations: 


XR  « kn^/(n2  + 1) 

XL  * kn^[  1 + 24(l-v)/(nL/rQ)2]/[n2  + 1 + 12/(nL/rQ)2]  (16) 

in  which  the  subscripts  "R"  and  "L"  of  XR  and  X^  identify  "Rayleigh" 
and  ’love",  respectively. 
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Discussion  of  Results 

The  geometric  and  physical  parameters  of  a family  of  cylindrical 
shells  selected  for  numerical  study  in  the  present  work  are  as  follows: 

L/rQ  = 5.33*  rQ/h  = 375 » v = 0.37.  These  values  were  chosen  according  to 
the  data  of  models  fabricated  for  a concurrent  experimental  study. 

Similar  to  the  procedure  adopted  in  Ref.  1,  the  first  phase  of 
computation  was  carried  out  to  determine  the  normal  modes  of  vibration  of  a 
circular  cylindrical  shell.  To  this  end,  modes  with  circumferential  wave 
number  0 < n < 29  were  obtained  for  the  first  four  axial  wave  numbers.  Here, 
the  symbol  m is  chosen  to  denote  such  a number  and  m has  the  assigned 
values  to  1 to  4.  In  addition,  a second  nominal  axial  wave  number  m is 
defined  to  me  m = m - 2.  Thus,  m = 1 and  2 represent  the  first  longitudi- 
nally symmetric  (Rayleigh)  modes  and  the  first  longitudinally  antisymmetric 
(Love)  modes,  respectively;  while  m=  3 and  4,  respectively,  commonly 
recognized  as  m = 1 and  2,  and  the  second  lowest  longitudinally  symmetric 
and  longitudinally  antisymmetric  modes.  It  should  be  noted  that  the  alter- 
nating alignment  of  symmetry  and  antisymmetry  of  modes  has  not  been  assumed. 
The  question  of  circumferential  symmetry,  of  course,  does  not  arise  for  the 
case  of  circular  cylinders. 

As  mentioned  in  the  preceding  section,  the  Rayleigh  modes  lie  in 
2 2 

zone  2 (&^  , dp  >0),  and  the  Love  modes  were  found  in  both  zone  2 and  3 
2 2 

(8j  < 0,  Sg  >0).  Corresponding  results  are  illustrated  in  Fig.  2,  which 

2 2 

displays  the  values  of  Bj  and  Bg  versus  the  circumferential  wave 

2 2 

number  n.  The  other  two  complex  conjugate  roots  B^  and  are  not 

2 2 

shown,  but  they  are  very  much  greater  than  Bj  and  in  magnitude. 

Thus,  Bj  and  Bg  characterize  the  long  wave  form  of  the  shell  deformation, 
which  combines  with  the  short  wave  form  characterized  by  Bj  and  8^. 


Although  an  extremely  high  degree  of  accuracy  is  required  for  the  eigen- 
value X in  order  to  accurately  determine  the  mode  shape,  the  trial  and 
error  procedure^  for  the  solution  of  the  frequency  equation  can  accomplish 

Y 

this  end  with  relative  ease.  In  all  cases  computed,  the  eigenvalue  X is 
determined  up  to  13  significant  figures.  With  such  accuracy  easily  attainable, 
the  closeness  of  the  first  two  eigenvalues,  due  to  the  narrowness  of  zone 
2,  does  not  cause  insurmountable  difficulties.  The  approximate  values  of 
X^  and  X^,  given  by  Eqs.(l6),  provide  valuable  guidance  to  the  determina- 
tion of  more  accurate  values  of  X.  For  example,  when  n = 9,  Eqs.(l6)  yield 
X^  « 3 .84058 . . .xl 0 ^ , and  X^  « 3<8654...xl0  ^ . The  corresponding  roots  for 
n = 9 and  m = 1 ,2  are  found  to  be  3 .8283 . . .xl 0-^  and  3 .8424 . . .xl 0-^ , 
respectively.  On  the  other  hand,  since  the  latter  two  values  of  X , close 
as  they  are,  yield  entirely  different  mode  shapes,  the  Rayleigh-Love  approx- 
imate frequencies,  i.e.,  X^  and  X^,  can  not  be  relied  upon  to  provide 
serviceable  mode  functions.  Since  the  roots  of  X determined  are  so  close 
to  one  another,  plots  such  as  that  partially  shown  in  Fig.  2 are  quite  val- 
uable in  the  process  of  identifying  m. 

Figures  3a  and  3b  show  the  mode  shapes  of  m = 1 and  m = 2, 
respectively,  for  some  selected  values  of  n , based  upon  solutions  of  Eqs .(12) 
to  (14).  They  demonstrate  that  the  deformation  of  the  former  is  longitudi- 
nally symmetric  and  that  of  the  latter  antisymmetric.  Such  a property  has 
not  been  assumed  a priori  in  the  analysis.  Moreover,  these  mode  functions  turn 
out  to  be  quite  different  from  the  Rayleigh-Love  approximation.  An  exam- 
ination of  the  numerical  amplitudes  of  the  functions  Un  and  Vn>  although 
smaller  than  that  of  the  corresponding  W^,  suggests  that  the  deformations 
of  the  first  two  families  are  not  truly  inextens iona 1 . 


Upon  the  determination  of  30  modal  functions  for  each  m = 1,2, 3, 4 
(15  even  and  15  odd;  i.e.,  N = 15),  they  are  substituted  into  the  variational 
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functional  5F,  and  the  required  numerical  integrations  with  respect  to 
x are  performed.  Lagrange  multipliers  are  introduced  to  incorporate  the 
2N  constraint  conditions  (Eqs . (9)  and  (11))  into  F,  and  are  subsequently 
eliminated.  After  all  these  straight-forward  steps  are  completed,  the  system 
of  equations  is  finally  expressible  in  the  form 

Pz  - XQz  = 0 (17)  .. 

of  which  the  dimensionality  is  reduced  to  N.  The  eigenvalues  and  the  cor- 
responding eigenvectors  are  determined  from  Eq .(17)  Results  are  shown  in 
Figs . 4 to  6 . 

In  Figs.4and  5 the  nond imens ional  natural  frequency  corre- 
sponding to  the  first  four  axial  wave  numbers  m = 1,2, 3, 4 are  shown  versus 
the  nominal  circumferential  wave  number  n.  The  resulting  variation  of 
deflection  along  the  circumference  of  a vibrating  oval  cylindrical  shell  can 
be  substantially  different  from  the  simple,  sinusoidal  pattern  of  an  equiv- 
alent circular  cyl indrical  shell.  Thus,  n cannot  be  determined  by  the 
count  of  nodes  along  the  circumference.  Rather,  the  determination  of  n is 
made  possible  by  performing  the  calculations  for  a series  of  values  of  the 
eccentricity  parameter  ? with  a small  increment  of  5 (chosen  to  be  0.1). 

' 

The  evaluation  of  the  mode  shape  can  then  be  traced  by  means  of  a plot  of 
the  frequency  versus  § exemplified  by  Fig.  6.  Tha  designation  of  a number 
for  of  a mode  n is  aided  by  the  continuous  variation  of  such  a plot. 

It  turns  out,  in  the  free-free  case,  that  the  values  of  n coincide  with 
the  index  of  the  largest  component  (in  absolute  value)  in  the  associated 
eigenvector  z ; such  is  not  always  the  case  for  the  supported  shells  studied 
in  Ref  . 1 . 

Figure  5 magnifies  the  region  of  low  frequency  and  small  n in 

I 

J 


which  the  contrast  between  the  spectra  of  the  circular  and  noncircular 
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cyl  inders  is  most  noticeable. 

It  is  significant  to  note  that,  for  any  mode  of  deformation  of 
an  oval  shell  represents  by  a point  in  Figs.  4 or  5»  a second  pattern  always 
exists.  More  precisely,  for  the  same  circumferential  wave  number  n,  a 
higher  frequency  can  always  be  found  depending  upon  the  status  of  symmetry 
about  the  principal  axes.  Similar  characteristics  were  observed  by  Culberson 
and  Boyd  in  their  investigation^  from  the  odd  modes  of  a freely  supported 
oval  shell,  and  by  the  present  authors'  from  both  the  odd  and  even  modes  of 
the  seven  other  supported  cases.  However,  unlike  what  was  observed  from  most 
of  the  supported  cases,  a pattern  seems  to  emerge  from  the  present  results 
in  that  a lower  frequency  is  associated  with  the  mode  shape  antisymmetric  with 
respect  to  the  major  axis.  In  terms  of  the  two  sets  of  modal  expansions 
( Eq s . (5 ) and  their  antisymmetric  counterparts)  and  the  sign  of  £ which  fixes 
the  origin  of  the  circumferential  coordinate  s , the  lower  mode  possesses 
the  following  specification: 

. (a)  Odd  modes:  The  modal  expansion  (lower  mode)  is  symmetric 

if  the  eccentricity  parameter  § is  chosen  to  be  negative;  and  it  is  anti- 
symmetric if  ? is  positive. 

(b)  Even  modes;  The  modal  expansion  (lower  mode)  is  antisymmetric 
regardless  of  the  sign  of  |. 

A visual  version  of  these  rules  is  given  in  Fig.  7-  In  Figs.  4 to  o,  only 
the  lower  frequencies  are  shown. 

Another  interesting  feature  of  the  free-free  shell  is  that  the 
dependency  of  the  natural  frequency  upon  the  eccentricity  parameter  is  not 
as  strong  as  in  the  supported  cases.  This  weaker  dependency  is  the  direct 
result  of  a weaker  coupling  among  the  circumferential  modes,  which  is  re- 
flected by  the  sharper  domination  by  one  component  in  any  eigenvector.  For 
this  reason,  only  curves  corresponding  to  the  extreme  values  of  the 
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eccentricity  parameter  |f|  = 0 and  1 are  shown  in  Figs.  4 and  5,  although 
computations  were  performed  for  all  intermediate  values  with  an  increment 
of  § equal  to  0.1.  Among  the  four  groups,  the  Rayleigh  modes  (m  = 1) 
exhibit  the  strongest  dependency  upon  I-  and  the  Love  modes  (m  = 2)  display 
the  weakest  dependence.  The  variation  of  the  natural  frequency  of  the 
Rayleigh  modes  (m  = 1)  versus  the  eccentricity  parameter  are  shown  in  Fig.6. 
Similar  curves  for  other  values  of  m,  while  not  being  shown  here,  exhibit 
both  increasing  and  decreasing  variations  of  the  frequency  as  !^|  increases, 
but  show  much  less  variation  with 

In  conclusion,  it  is  observed  that  the  rigorous  enforcement  of 
boundary  conditions  in  the  dynamic  analysis  of  cylindrical  shells  of  finite 
length  is  essential.  Moreover,  the  present  method  of  analysis  succeeded  in 
transforming  a two-dimensional  (spatial)  problem  into  a one-dimensional  one. 
It  seems  reasonable  to  predict  that  the  convenient  features  of  the  present 
method  are  equally  applicable  to  problems  involving  forced  dynamic  response 
of  noncircular  cylindrical  shells.  Finally,  a remark  of  caution  regarding 
the  completeness  of  the  modal  functions  in  the  functional  basis  must  be 
made.  For  any  given  value  of  m,  granting  truncation,  the  complete  set  of 
modal  functions  must  be  included  in  the  modal  expansion.  A missing  link 
would  severely  damage  the  accuracy  of  the  neighboring  modes  as  a result  of 
coupl i ng . 
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FIG.  2 CHARACTERISTIC  WAVE  NUMBERS  (LONG  WAVES)  OF  THE  FIRST  AND 

SECOND  MODES  OF  A CIRCULAR  CYLINDRICAL  SHELL  WITH  FREE  ENDS. 
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NATURAL  FREQUENCY  OF  FREE-FREE  OVAL  SHELLS  VS.  CIRCUMFERENTIAL 
WAVE  NUMBER  (L/r  = 5.33,  r _/h  = 375,  v = 0.37)  . 


